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On  a Particular  Life  Testing  Problem 
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A.  John  Petkau* 


1.  INTRODUCTION 

Consider  the  usual  life  testing  situation  in  which  observations 
become  available  in  an  ordered  manner.  In  this  situation,  it  is  natural 
to  consider  the  possibility  of  terminating  the  life  test  at  an  early 
stage;  say,  after  the  first  r out  of  the  total  of  n observations 
have  been  recorded.  Any  such  procedure  has  the  advantage  that  it  may 
lead  to  a decision  in  a shorter  period  of  time  and  with  fewer  observations 
than  a procedure  which  involves  observing  what  happens  to  all  the  items 
being  tested. 

The  special  case  in  which  the  underlying  lifetime  distribution 

is  assumed  to  be  exponential  with  unknown  location  parameter  £ and 

unknown  scale  parameter  0 has  recently  been  considered  by  Perng  (1977). 

He  argues  that  a problem  which  will  sometimes  arise  in  practice  is  that 

of  testing  the  simple  null  hypothesis  H.  ; 0 e versus  the  composite 

u u 

alternative  H,  : 6 c £1  - ft.  where  0 * (3,0),  il  « {0  | CKoso,,  d-*3  ' 
10  0 0 

and  0,  » {0  I c-o,,  3-3,}  . Without  loss  of  generality,  we  will  assume 
0 0 u 

from  this  point  on  that  3^-0  and  a - 1 . 

Perng  argues  (Perng,  1977,  p. 1401)  that  'the  likelihood  ratio 
method  does  not  give  us  any  "neat"  test  for  this  problem,  since  the 
likelihood  ratio  test  statistic  is  rather  complicated,  and  the  distribution 

of  the  likelihood  ratio  test  statistic  under  H is  not  known'  and 

0 

consequently  13  motivated  to  look  for  other  approaches  to  the  problem. 

He  proposes  a test  procedure  based  on  Fisher's  method  of  combining 

independent  test  statistics  (see,  for  example,  Fisher  1950,  op.  19-101) 
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sense  of  Bahadur  efficiency. 

In  Section  2 of  this  paper  Perng's  test  procedure  is  described 
in  detail.  Critical  values  for  this  test  are  available  and  the  test  • 

appears  to  be  straightforward  to  carry  out.  Unfortunately,  as  is  pointed 
out  in  Section  2,  difficulties  arise  in  attempting  to  carry  out  the  test 
in  practice.  An  alternate  test  procedure  is  proposed  which  is  truly 
straightforward  to  carry  out.  Petkau  (1978)  has  already  shown  that  this 
alternate  test  is  asymptotically  suboptimal  in  the  sense  of  Bahadur 
efficiency.  Since  neither  of  these  alternate  test  procedures  is  entirely 
satisfactory,  we  return  to  consideration  of  the  likelihood  ratio  method. 

The  likelihood  ratio  test  statistic  is  determined  and  the  appropriate 
critical  values  for  carring  out  the  test  are  provided.  In  addition,  the 
asymptotic  null  distribution  of  this  test  statistic  is  determined.  The 
critical  regions  and  the  power  functions  of  these  three  procedures  are 
briefly  compared  in  Section  3. 


i 


2.  TEST  PROCEDURES 


Suppose  X .X  .-..X  are  the  order  statistics  of  a random 
sample  of  size  n from  the  exponential  distribution  with  probability 
density  function 


Vz) 


1 - (z-6) /a 

— e 

a 


for  z > 8 , 


where  9 * (8, a)  . We  consider  test  procedures  based  on  the  first 
r(l<r£n)  order  statistics  for  testing 


: 9 £ versus  H.  : 9 E ft  - ft n , 

0 0 1 o 


where  fl  * {(8,0)  | 8£0,  0<asl}  and  (2^  = {(0,1)}  . 

Define  the  statistics  V and  U as  follows 

n n 

r 

V = 2 f y (X.-X  ) + (n-r)  (X  -X.)]  , 
n i r r 1 


U = 2nX,  . 
n 1 


(2.1) 


In  the  general  case,  V /o  and  2n(X1~S)/a  are  independently  distributed 
2 2 

as  X(2r  2)  411(1  X(2)  respectively  and,  thus,  when  HQ  is  true,  the 

same  is  true  of  V and  U . Small  observed  values  of  V and  large 
n n n 

observed  values  of  U constitute  evidence  against  H . 

n 0 

Perng  (1977,  p.  1402)  defines  the  statistic 


where  T 

n 


Q = -2  loc  (U  -u)  -2  log  P„  (T  >t) 

n ' n0  n Ho  n 

V-'1'  and  u , t are  the  observed  values  of  V 
n n 


(2.2) 


and  T 

n 
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respectively  and  proposes  the  test  procedure 


Reject  Hq  if  i c^ 


Since  Q is  distributed  as  x,,,  under  , the  critical  value 

n (4)  0 

c_  * x'l  , can  he  determined  from  a chi-square  table.  However,  in 
Q 4,1-a 

order  to  carry  out  this  test  procedure,  one  must  evaluate  both  of  the 

probabilities  appearing  in  the  definition  of  the  statistic  . 

P„  (U  5u)  can  be  easily  determined  since,  under  , U /2  is  exponentially 
n un 

distributed  with  mean  one.  On  the  other  hand,  since  chi-square  tables 

contain  only  a few,  and  typically  very  few,  selected  precentiles,  the 

probability  P (T  2t)  cannot,  in  general,  be  accurately  determined. 

H0  n 

This  difficulty  is  familiar  to  even  the  most  casual  user  of  statistical 
methods  who  may  at  some  time  have  attempted  to  determine  an  observed 
P-value  from  a chi-square  table.  Approximations,  such  as  the  Wilson- 
Hilferty  approximation,  might  serve  to  alleviate  this  difficulty  some- 
what (for  interesting  related  work,  see  Hoaglin  (1977)),  however,  it  is 
clear  that  while  Perng's  test  is  easy  to  carry  out  in  theory,  difficulties 


arise  in  attempting  to  carry  out  the  test  in  practice. 

Since  V and  U (see  (2.1))  are,  under  H.  , independently 
n n 0 

2 2 

distributed  as  x,-.  and  X,-,%  , . ,,  , 

(2r-2)  A(2)  respectively,  and  since  small  values 


of  V and  large  values  of  U are  evidence  against  , an  alternate 

n n 0 

test  which  immediately  suggests  itself  is  to  define  the  statistic 


F = (r-l)U  /V 
n n n 


and  to  use  the  procedure 


...  .<r 
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Reject  H if  F > c 
J 0 n F 


Under  , F has  a F distribution  with  2 and  2r  - 2 degrees  of 
u n 

freedom  and  the  critical  value  c_  = F_  _ _ , can  be  determined  from 

F 2,2r-2,l-a 

a table  of  the  F distribution.  This  alternate  test  is  truly  straight 
forward  to  carry  out.  Unfortunately,  at  least  for  large  sample  sizes 
this  test  will  perform  poorly  in  comparison  to  Perng's  test  since,  as 
has  been  shown  by  Petkau  (1978) , this  alternate  test  is  asymptotically 
suboptimal . 

Since  neither  of  the  test  procedures  considered  above  is 

entirely  satisfactory,  we  now  consider  the  likelihood  ratio  test  procedure. 

★ 

In  what  follows,  denotes  the  vector  of  the  first  r order  statistics. 

The  marginal  joint  density  of  the  first  r order  statistics  (Epstein 
and  Tsao,  1953)  is  given  by 


gQ (x  ) = 
9 r 


(n-r) !o‘ 


exp{-[  y (x. -8)  + (n-r)  (x  -8) ] /a } 

r 


for  x > 3 , 


where  0 = (8, a)  . Defining 


vv 


sup  g. (x  ) 
9 r 

★ 

sup  g (x  ) 
9efl  9 r 


it  is  easy  to  verify  that  this  likelihood  ratio  statistic  can  be  written 
(see(2.1))  as 


■ 


i. 


- 6 


X (X  ) * exp  (-(J  /2) 
nr  n 


exp (-U  /2) { (V  /2r)exo[l-(V  /2r)]}r 
n n n 


for 

V > 
n 

2r  , 

(2.6) 

for 

V s 
n 

2r  . 

The  usual  likelihood  ratio  test  procedure  is 


Reject  H if  X £ k 
J 0 n L 


(2.7) 


Perhaps  a more  convenient  form  is  given  by 


- 2 log  X (X  ) * U 
nr  n 


for  V > 2r  , 
n 


(2.8) 


U - 2r{ 1 - (V  /2r)  + log(v  /2r) } 
n n n 


for  V £ 2r  , 
n 


with  the  corresponding  procedure 


Reject  H if  -2  log  \ > -2  log  k = c 

0 n L L 


(2.9) 


It  should  be  noted  that  this  test  statistic  is  easy  to  evaluate  and  it 

reamins  only  to  evaluate  the  critical  value  k , or  equivalently,  c 

I*  L 

In  what  follows  p (x)  , f (x)  denote  the  densitv  function  and  distribution 
mm 

function  respectively  of  a random  variable  having  a gamma  distribution 
with  parameters  m and  1 , that  is. 


P ( x ) 
m 


T (m) 


m-1  -x 
x e 


for  x > 0 , 


V (x) 
m 


ip  (y)dv  . 
m 


0 


I 


* ' 
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We  also  recall  that  if  2 is  distributed  as  x,,  » then  z/2  is 


(2m) 


distributed  as  a gamma  with  parameters  m and  1 . The  following  theorem 
evaluates  the  critical  value  k appearing  in  (2.7). 

L 


Lemma  2.1.  The  appropriate  critical  value  k^  for  a size  a (0<a<l) 


test  is  given  by 


1-x 

kL  “ [V  ]r 


(2.10) 


where  is  the  unique  solution  between  0 and  1 of 


l-x,r 


[xe"  ] {l-'?r_1(r)  + (l/x-l)ri()r_1(r)  } + <?  (rx)  = a . (2.11) 


Proof:  For  0 < t < 1 , we  have 


PH  (WSt)  " PH  ('log  W " ' log  ' 
0 0 


j P (-log  X (x  ) > - log  t | V /2r  = x)dPu  (V  /2r  s x)  , 
J n r n H0  n 


which,  by  (2.11)  ana  the  independence  of  U and  V under  H , can  be 

n n 0 


for  x between  0 and  1 , there  is  a unique  solution,  xQ  , between 
0 and  1 of  the  equation  h(x)  =0  which  can  be  written  as 


r l-x,r 
t = [xe  ] . 


(2.12) 


Then 


P (X  (X)St) 
H0  n r 


IdP  (V  /2r^x)  + P (U  /25-log  t+r[l-x+log  x])dP  (V  /2rix) 

J fi  n J H0  n H0  n 

0 X0 


+ P„(U  /25-log  t)  dPti  (V  /2r<x)  , 
H„  n I H.  n 

0 1 0 


= f , (rx„)  + t re 
r-1  0 


-r 


x rerX,^r  ^(rx)dx 


+ t[l  - 4'r_ L ( r ) ] , 


which  can  be  written  as 


p„  (X  (X  ) St)  = V , (rx  ) + t{l-f  , (r)  + r (1/x  -1)  ip  _ (r)  }.  (2.13) 

HQ  n r r-1  0 r-1  0 r-1 


This  gives  the  null  distribution  function  of  the  likelihood  ratio  statistic 
which  leads  directly  to  the  desired  result. 

It  should  be  noted  that  the  critical  values  k , or  equivalently 


c = -2  log  k , are  functions  of  a and  r and  do  not  depend  upon  n . 

Li  Li 


f 


9 


Although  equations  (2.10)  - (2.11)  appear  rather  complicated,  the 
availability  of  standard  numerical  routines  whic  evaluate  the  incomplete 
gamma  integral  permits  the  easy  evaluation  of  the  critical  values  by 
means  of  an  iterative  procedure.  Indeed,  if  desired,  the  null  distribution 
of  the  likelihood  ratio  statistic  could  be  examined  in  detail.  Critical 
values  c^  for  selected  values  of  a and  r are  provided  in  Table  1. 

Table  1 about  here 


It  is  of  interest  to  examine  the  asymptotic  (r-**)  null 
distribution  of  minus  two  times  the  loo  of  the  likelihood  ratio  test 
statistic.  A cursory  glance  at  Table  1 reveals  that  convergence  to  this 
limiting  distribution,  whatever  it  be,  is  quite  slow.  Naively  one  might 
expect  this  limiting  distribution  to  be  chi-square,  as  is  customary. 

That  this  is  not  quite  the  case  is  established  in  the  following  theorem. 


Theorem  2.1.  The  asvmototic  (r—»)  null  distribution  of  -2  log  \ (X  ) 
nr 

2 2 

is  a mixture  of  and  » each  with  probability  1/2  . 


Proof:  Replacing  t by  e 


-t/2  in  (2.12'  and  (2.131,  we  have 


P (-2  log  \ (X  )it)  = Y (rx  ) + e t'  (r)+r(l/x  -IK'  , (r) '•  (2.1 


r-1  0 


r-1 


0 r-1 


for  t > 0 , where  x ^ is  the  unique  solution  between  0 and  1 of  the 
equation 


-t/2  r i-x.r 
e » .x  e 


(2.1? 


10 


In  order  to  evaluate  the  limit  as  r ■»  » in  (2.14),  we  examine  the 
behavior  of  x^  as  r ■*  20  . Writing  (2.15)  in  the  form 


1 - x + log  x = - t/2r  , 


it  is  clear  that  as  r •*  30  , xQ  -*■  1 . Simple  expansions  yield  that 
as  r •*  30  , 


x - 1 - A + A2/3  ♦ A3/12  + 


(2.16) 


1/2 

where  A = (t/r)  . In  the  following  denotes  a gamma  random 

variable  with  parameters  r - 1 and  1 and  2 = [w^  - (r-l)]/(r-l) 

It  is  well  known  that  Z has  an  asvmDtotic  distribution  which  is 

r 

standard  normal.  Examining  the  terms  in  (2.14),  we  have  as  r -*•  00  , 


1/2 


-1/2 

, (r)  = P(W  sr)  = P (2  <(r-l)  ) - i(0)  = 1/2  , 

r-1  r r 


l/o  l/o 

4'r-l(rx0)  = P(wr-rx0)  = p(2r  " [r (Xq-D  + l]/(r-l)  ) - *(-t  ' ) , 


where  in  the  latter  we  have  made  use  of  (2.16).  Further,  use  of  (2.16) 
together  with  Stirling's  formula  yields  that  as  r -*  °>  , 


r ( 1/x  -1)  (r) 

0 r-1 


Collecting  these  results,  we  obtain 


1/2  . , _ . 1/2 
t /(2tt) 


lim  ? ( -2  log  \ (X  ) 2 t)  * $(-t  ) 

r— ” ‘ 0 


1/2,  1 -t/2 


2 a 


1/2  -t/2 
t e 


1 -t/2 

2 6 


. 1/2  -t/2 

t e 


/2* 


which  is  the  desired  result. 

It  is  easy  to  verify  that  the  (l-oi)th  quantile  of  this 
asymptotic  distribution  is  given  by  c , say,  where  c/2  is  the 
solution  of  the  equation 


f3/2(x)  - e 


-x 


l-2oi 


(2.17) 


These  critical  values  are  displayed  for  selected  values  of  a in  the 

last  line  of  Table  1.  It  should  be  noted  that  since  the  convergence 

of  the  critical  values  c to  the  asymptotic  critical  values  c is 

L 

quite  slow,  this  asymptotic  distribution  is  not  particularly  useful  for 
practical  purposes.  Indeed  the  magnitude  of  the  effect  of  using  the 
asymptotic  critical  values  as  approximations  to  the  actual  critical  values 
can  be  roughly  estimated  for  any  fixed  r from  Table  1.  If  desired, 
this  effect  could  be  examined  in  detail  by  evaluating  (2.12)  - (2.13) 
with  exp (-c/2)  replacing  t . 

The  availability  of  critical  values  for  the  likelihood  ratio 
test  makes  this  procedure  easy  to  carry  out.  None  of  the  practical 
difficulties  which  arise  in  carrying  out  Perng's  test  are  encountered. 

From  this  point  of  view,  both  the  F test  and  the  likelihood  ratio 
test  are  preferable  to  Perng's  test.  On  the  other  hand,  Perng's  test 
is  known  to  be  asymptotically  optimal  in  the  sense  of  Bahadur  efficiency 
while  the  F test  is  asymptotically  suboptimal.  Further,  Bahadur  and 


12 


Raghavachari  (1972)  have  shown  that  any  test  procedure  based  on  a likelihood 
ratio  test  statistic  is,  under  regularity  conditions,  asymptotically 
optimal.  Although  these  asymptotic  considerations  rule  out  the  F test 
as  a viable  competitor  to  either  Perng's  test  or  the  likelihood  ratio 
test  for  all  large  enough  sample  sizes,  no  preference  is  indicated  between 
these  latter  two  tests.  Moreover,  the  implication  of  these  asymptotic 
considerations  for  the  fixed  sample  size  case  is  unclear.  In  the  next 
section  we  consider  the  fixed  sample  size  properties  of  these  three  test 
procedures . 


13 


3.  COMPARISON  OF  TEST  PROCEDURES 


Since  critical  values  for  the  likelihood  ratio  teat  are  now 


available,  this  procedure  would  seem  to  be  preferred  to  Perng's  teat 
in  view  of  the  practical  difficulties  encountered  in  carrying  out.  the 
latter.  On  the  other  hand,  since  Ferng's  test  is  asymptotically  optimal, 
it  is  of  interest  to  compare  the  behavior  of  Perng's  test  with  that  of 
the  likelihood  ratio  test  in  fixed  sample  sizes.  Further,  it  is  of 
interest  to  determine  whether  the  asymptotic  inferiority  of  the  F test 
is  also  evident  in  the  fixed  sample  size  case.  In  this  section  we 
examine  and  compare  the  critical  regions  and  the  power  functions  of  the 
three  teat  procedures. 

Recalling  (-.1),  we  now  define 


V - U /2  , 

n n 

Z - V /lx  . 
n n 


From  the  distribution  theory  of  Section  l,  we  have 


d 

V » nt>  + o • w 
n 


„ d o 
z • — • w 
n r r-1 


where  w and  w ^ are  independent  random  variables  having  (unit) 
exponential  and  gamma  (with  parameters  r-1  and  1)  distributions 
respectively.  While  the  distribution  of  Y is  the  same  for  all  n 
(except  for  the  shifting  of  the  origin  to  the  point  ni)  , the  distribution 


of  Z becomes  degenerate  at  o as  r . We  will  now  describe 

n 


the  critical  regions  of  the  test  procedures  of  Section  2 in  the  IY  ,Z  ) 

n n 

plane.  Note  that  Y and  2 are  both  positive  random  variables. 

n n 

From  (2.8)— (2.9) r the  critical  region  for  the  likelihood  ratio 
test  is  given  by  the  union  of  the  two  sets  A and  B defined  by 

A » {z>l}  n {y  i c/2)  , 

Li 

B - {zsl}  n (y  - r[l-z  + log  z]  i c /2 } . 

L 

Since  1 - z +•  log  z S 0 for  z S 1 , the  set  B contains  all  (y,z) 

pairs  for  which  z S 1 and  y z c /2  . Thus,  all  pairs  for  which 

Li 

y ^ c/2  and  no  pairs  for  which  both  v < c /2  and  z > 1 are  contained 
Li  L 

in  the  critical  region.  For  the  remaining  pairs,  namely  those  with 
z s 1 and  y < c /2  , only  those  points  below  (or  on)  the  curve  z = z(y) 

Lj 

defined  by 

y - r[l  - z + log  z]  = cL/2  (3.2) 

are  in  the  critical  region.  Examining  this  curve,  we  have  z(c  /2)  = 1 

L 

and  z(0)  = z (c  ,r)  , say,  is  the  unique  number  between  0 and  1 
L L 

which  is  the  solution  of 

r[l  - z + log  z]  + c /2  - 0 . 

L 

Further,  ^ = s/[r(l-z) ] which  is  positive  for  z < 1 and  approaches 

dy 

* as  z •*■  1 , that  is,  as  y -*■  c /2  . Thus  the  curve  z(y>  is  monotone 
increasing  for  0 s y < c /2  and  fits  smoothly  with  the  vertical  line 
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y * c /2  at  the  point  (c  /2,  1).  Indeed,  a simple  expansion  yields 

Li  Li 

that  for  z(y)  defined  by  (3.2), 

z(y)  ~ 1 - (2c/r) 1 


when  y ■ c /2  - e , where  € is  a small  positive  quantity. 

Li 

We  note  that,  for  fixed  r , z (c  ,r)  decreases  as  c 

Li  Li  Li 

increases,  that  is,  as  the  size  of  the  test  decreases.  Further,  for 
fixed  a , c ♦ c defined  by  (2.17)  as  r -*  * . Thus,  c./r  •*  0 

Li  Li 

as  r * • and  consequently  z (c  , r)  -*■  1 as  r * 00  . As  r increases, 

Li  Li 

the  critical  region  of  the  likelihood  ratio  test  approaches  the  set  of 
all  points  in  the  first  quadrant  not  satisfying  both  v i c/2  and  z i 1 . 

Remark:  Given  this  simple  limiting  form  of  the  critical  region,  at  first 
glance  it  may  seem  attractive  to  use  this  limiting  region  as  an  approximation 
to  the  critical  region  for  fixed  r . That  this  should  not  be  done  is 
clear  upon  evaluation  of  the  size  of  the  resulting  test 


P (reject)  ■ 1 - P (Y  S c/2) »P  (2  2 1) 

H0  H0  n Ho  n 

* 1 - (1  - e‘c  2) • (1  - f , (r) ) 

r-1 


Now  V , (r)  is  monotone  decreasing  in  r (for  r 2 2)  and  V , (r)  -*■  1 2 
r-1  r-1 

as  r ■*  » . Thus,  for  the  test  which  corresponds  to  this  critical  region 
we  have 


1 - (1-e  ) ( 1-e  ) 2 P (reject)  2 1 - (1-e 

**  a 


a-0  /2  2 1/2 


Turning  to  consideration  of  Ferng's  test,  we  see  from  (2.2) 


(2.3)  that  the  critical  region  is  the  set  of  all  iy,z)  pairs  satisfying 


y - log  ¥ (r*z)  5 c /2  . 


Since  - log  ¥^  ^(r*z)  > 0 for  z > 0 , all  (y,z)  pairs  for  which 
y 2 Cg/2  are  contained  in  the  critical  region.  Examining  the  curve 
z - z(y)  defined  by  equality  in  (3.3),  we  have  z(c^/2)  « « and 
z(0)  * z^(c^,r)  , say,  is  the  unique  positive  number  satisfying 


- log  ¥ (r*z)  = c /2  . 

r-1  Q 


Recall  that  c = v"  , • Since  ¥ , (r)  is  monotone  decreasing 

Q r,l-a  r-1 

and  approaches  1/2  as  r -*•  » , z^(c^,r)  < 1 provided  that 

c^/2  - log  2 > 0 which  will  certainly  be  the  case  for  any  a of  interest. 


Further , 


¥ (rz) 
dz (y)  r-1  r 


dy  ry  , (rz) 

r-1 


which  is  positive  and  approaches  00  as 


z -*■  » , that  is , as  y -*■  c^/2  • Thus,  the  critical  region  consists  of 
those  points  below  (or  on)  the  curve  z(y)  which  is  monotone  increasing 


and  increases  without  bound  as  v ■*  c /2  . 

Q 


We  note  that,  for  fixed  r , z (c  ,r)  decreases  as  c increases 

Q Q Q 

that  is,  as  the  size  of  the  test  decreases.  Further,  it  is  clear  that, 
for  fixed  at  , z^(c  , r)  1 as  r -*■  » . For  purposes  of  comparing  the 

X X 

critical  region  of  Perng's  test  to  that  of  the  likelihood  ratio  test,  it 

★ 

is  of  interest  to  determine  the  point  (y  , 1)  which  is  on  the  curve  defined 
by  (3.3).  We  have 


y (c  ,r)  = log  ¥ ,(r)  + c'2 

Q r-1  Q 


We  assume  as  before  that  c/2  - log  2 > 0 which  insures  that  the  right 

V 

¥ . (r)  - 1/2 

r-1 


hand  side  is  positive.  As  r increases, 


and  we  have 
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these  procedures/  it  would  clearly  be  of  interest  to  obtain  explicit 
expressions  for  the  power  functions  of  these  procedures.  These  power 
functions  will  in  each  case  depend  upon  the  number  of  observations,  n , 
and  the  number  of  order  statistics  being  observed,  r , in  addition  to 
the  parameter  pair  (S,o)  and  the  size  of  the  test.  The  straightforward 
calculations  leading  to  these  power  functions  are  not  displayed;  the 
results  are  merely  summarized  in  the  following  lemma. 


Lemma  3.1.  (a)  Likelihood  ratio  test: 

Set  A = (c  /2  - n6)/0  . The  power  function,  P , 

LR 

is  given  by 


P =1 
LR 


for  A s 0 


-A 


e [1  - ,i'r_1(r/o)  ] + fr_1(rh/a)  + 


(3.5) 


( r—  1 ) 


[r(l-a)/0+l] 


e‘\r(r/a)[h'trll-a)/a+1]  - l] 


for  A > 0 , 

where  h is  the  unique  solution  between  0 and  1 of 

r[l  - h + log  h]  * - a A . (3.6) 


(b)  Preng's  test: 

Set  y = (c  /2  - n8)/a  . The  cower  function,  ° 

Q * PR 

is  given  by 


Set  £ 3 a^'/^r  ^ and  n - nB/j  . The  oower  function, 
Pp  , is  given  by 


poorly  relative  to  the  other  tests  for  alternatives  where  n 


nB/o 


! 


I 


vTTu^;..  I 
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is  small  and  a is  also  small.  The  rather  complicated  form  of  the 
expressions  appearing  in  Lemma  3.1  makes  comparison  among  these  power 
functions  difficult  other  than  by  means  of  evaluation  for  specific 
choices  of  a , r , a and  nB  . 

Equations  (3.6)  and  (3.8)  have  already  been  encountered  and 
solving  these  equations  presents  no  particular  difficulty.  The  power 
functions  of  each  of  these  test  procedures  can  then  be  evaluated  with 
the  help  of  a standard  numerical  routine  which  evaluates 
the  incomplete  gamma  integral  provided  that  the  integral  appearing  in 
(3.7)  can  be  evaluated.  In  this  connection  note  that 


-1 

a 


/■« 


<i'r-l(w/a) 

[fr-1(w)]1/a 


dw  S [1  - ¥ (A/O) ]/[Wr_  (A) ]1/0  . 


For  fixed  r and  a , this  bound  is  monotonically  decreasing  in  A and 
consequently,  for  large  enough  A , this  integral  is  negligible.  Thus, 
the  integral  appearing  in  (3.7)  can  easily  be  evaluated  by  any  of  the 
standard  routines  which  numerically  integrate  a prescribed  function  over 
a finite  range. 

Tables  of  the  power  functions  (3.5),  (3.7)  and  (3.9)  have  been 
evaluated  for  various  choices  of  r and  a for  grids  of  (n 3, a)  pairs. 
It  is  rather  difficult  to  summarize  these  results  concisely  but  it  does 
appear  that  the  likelihood  ratio  test  is  more  sensitive  to  departures  in 
location  than  is  Perng’s  test.  More  precisely,  the  only  region  of  the 
parameter  space  in  which  the  power  functions  of  Perng’s  test  and  the 
likelihood  ratio  test  seem  to  differ  substantially  is  the  region  of  a 


close  to  1 and  nB 


reasonably  large  where  the  likelihood  ratio  test 


■ 


aa 
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has  a distinct  advantage.  In  fairness,  we  note  that  Perng's  test  appears 
to  be  marginally  more  sensitive  to  departures  in  scale.  As  already 
noted,  for  fixed  r and  a , the  power  function  of  the  F test  depends 
only  upon  n = nS/a  and  consequently  behaves  quite  differently  from 
the  power  functions  of  the  other  two  tests.  It  is  interesting  to  note  that  th 
power  function  of  the  F test  tends  to  dominate  the  power  functions  of 
both  of  the  other  tests  for  moderate  values  of  nS  and  the  difference 


can  be  substantial.  On  the  other  hand,  as  was  anticipated,  for  the  region 
in  which  both  nB  and  o are  small,  the  F test  performs  poorly  relative 
to  the  other  two  tests.  The  powers  for  the  case  a = 0.05  and  r = 10, 

20  and  40  are  displayed  for  selected  values  of  nS  and  a in  Table  2. 


Table  2 about  here 


4 . CONCLUDING  REMARKS 


We  have  considered  three  different  test  procedures  for  the 
particular  life  testing  problem  considered  by  Perng  (1977).  Critical 
values  which  make  the  likelihood  ratio  test  easy  to  carry  out  have 
been  provided.  Comparison  of  the  power  functions  of  Perng 's  test  and 
the  likelihood  ratio  test  reveals  that  although  these  power  functions 
are  quite  similar  over  most  of  the  parameter  space,  the  likelihood  ratio 
test  has  the  advantage  in  the  region  of  the  parameter  space  where  there 
is  a substantial  difference  in  these  power  functions.  Thus,  the  likeli- 
hood ratio  test  should  be  preferred  to  Perng' s test.  This  preference 
is  strengthened  in  view  of  the  practical  difficulties  associated  with 
carrying  out  Perng 's  test.  Examination  of  the  power  function  of  the  F 
test  reveals  that  in  the  fixed  sample  size  case  it  is  not  possible  to 
eliminate  the  F test  as  a viable  competitor  to  the  other  two  tests 
since  there  is  a region  of  the  parameter  space  in  which  the  F test 
does  substantially  better  than  either  the  likelihood  ratio  test  or 


Perng' s test. 
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Table  1.  Distribution  of  -2  log  X. 
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Tabl«  2.  (cont'd) 
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Table  2.  (cont’d) 
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For  a particular  hypothesis  testing  problem  involving  the 
location  and  scale  parameters  of  an  exponential  distribution,  Perng 
(1977)  proposed  a test  procedure  based  on  the  first  r out  of 
n observed  failure  times.  In  this  paper,  the  likelihood  ratio  Test 
for  the  same  problem  is  considered.  Critical  values  are  provided  and 
the  asymptotic  null  distribution  is  determined.  An  alternate  test 
based  on  an  F statistic  is  also  proposed  and  the  critical  regions 
and  power  functions  of  all  three  procedures  are  examined. 
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